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Abstract 

There have been many theoretical models constructed which aim to explain the neutrino masses and 
mixing patterns. While many of the models will be eliminated once more accurate determinations of the 
mixing parameters, especially sin^20i3, are obtained, charged lepton flavor violation experiments are able 
to differentiate even further among the models. In this paper, we investigate various rare lepton flavor 
violation processes, such as ii — > + 7 and ii — e conversion, in five predictive supersymmetric (SUSY) 
5'0( 10) models and their allowed soft-SUSY breaking parameter space in the constrained minimal SUSY 
standard model. Utilizing the Wilkinson Microwave Anisotropy Probe dark matter constraints, we obtain 
lower bounds on the branching ratios of these rare processes and find that at least three of the five models 
we consider give rise to predictions for /i — > e + 7 that will be tested by the MEG Collaboration at PSI. 
In addition, the next generation /i — e conversion experiment has sensitivity to the predictions of all five 
models, making it an even more robust way to test these models. While generic studies have emphasized 
the dependence of the branching ratios of these rare processes on the reactor neutrino angle, ^13, and the 
mass of the heaviest right-handed neutrino, Afa, we find very massive M3 is more significant than large 613 
in leading to branching ratios near to the present upper limits. 
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I. INTRODUCTION 

The recent compilation of oscillation data from the atmospheric [l| , reactor , and long baseline 
Q] neutrino experiments has provided solid evidence that neutrinos have small but non-zero masses. 
A global fit to current data gives the following 2a limits for the mixing parameters 

sin^ 6*12 = (0.28 - 0.37), Am|i = (7.3 - 8.1) x 10"^ eV^ (1) 
sin^ (923 = 0.38 - 0.63, Amj^ = (2.1 - 2.7) x 10"^ eV^ (2) 
sin^ 6i3 < 0.033 . (3) 

Since then, the measurements of neutrino oscillation parameters have entered a precision era. On 
the other hand, as no information exists for the value of ^13, the Dirac or Majorana nature of 
the neutrinos, the Dirac and/or Majorana CP phases, and the neutrino mass hierarchy, there are 
discoveries that are still yet to come. 

In the Standard Model (SM), due to the lack of right-handed neutrinos and the conservation of 
lepton numbers, neutrinos are massless. To generate non-zero neutrino masses thus calls for physics 
beyond the SM. There have been many theoretical ideas proposed with an attempt to accommodate 
the experimentally observed small neutrino masses and the larger mixing angles among them. In 
Ref. [S], we have surveyed 63 models in the literature that are still viable candidates and have 
reasonably well-defined predictions for ^13 . We found that the predictions for sin^ 2^13 of half of 
the models cover the range from 0.015 to the present upper bound of 0.13. Consequently, half of 
the models can be eliminated in the next generation of reactor experiments. 

One of the implications of the observation of neutrino oscillation is the possibility of measurable 
branching ratio for charged lepton flavor-violating (LFV) decays. While not the case in the SM, 
predictions of the supersymmetric (SUSY) Grand Unified Theories (GUT) for these rare decays are 
much enhanced, as these processes are suppressed by the SUSY scale, rather than the Plank scale 
Q]. Furthermore, as different models obtain large neutrino mixing angles through different mech- 
anisms, their predictions for the LFV charged lepton decays can be very distinct. Consequently, 
LFV charged lepton decays may provide a way to distinguish different SUSY GUT models. 

Among the models aiming to explain the neutrino masses and mixing, a particularly promising 
class are those based on (SUSY) SO(IO); for recent reviews of SO(IO) models, see In this 
paper, we investigate the predictions for various LFV charged lepton decays as well as muon- 
electron conversion in five of the SUSY SO(IO) models, assuming constrained minimal SUSY 
standard model (CMSSM) boundary conditions where only five soft-SUSY breaking parameters 
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are present. Furthermore, we impose the Wilkinson Microwave Anisotropy Probe (WMAP) dark 
matter constraints in the neutrahno, stau and stop coannihilation regions. Specifically, we present 
the allowed values for these soft-SUSY parameters for various branching ratios of these rare LFV 
processes. In addition, the lower bounds on the predictions for these rare processes in the five 
different SUSY SO(IO) models are given. We find that the predictions in these models are very 
distinct. We note the crucial role of the WMAP constraints in deducing the lower bounds on the 
predictions. 

Many authors have previously studied the branching ratio predictions for charged LFV decays 
in the SUSY GUT framework. Rather than study specific models, they generally adopt a generic 
approach and assume a nearly diagonal Cabibbo-Kobayashi-Maskawa (CKM)-like or bimaximal 
Pontecorvo-Maki-Nakagawa-Sakata (PMNS)-like mixing matrix to diagonalize the Yukawa neutrino 
matrix ^|. Following the procedure of Casas and Ibarra to invert the seesaw formula, they have 
carried out Monte Carlo studies by scanning the unknown right-handed neutrino mass spectrum and 
the angles and phases of the inversion matrix in order to present scatter plots of the rare branching 
ratios. A few exceptions to this procedure can be found in Ref. [l^. Here we are interested in 
specific models in order to determine the ability of the LFV experiments to differentiate among 
and rule out some models. The models chosen are highly predictive and illustrate a wide variety of 
textures for the charged lepton and neutrino mass matrices, and predictions for the right-handed 
neutrino mass spectrum and the reactor neutrino mixing angle, ^13. 

In Sec. II we describe the five representative SUSY SO(IO) models that we have analyzed. 
In Sec. Ill we review LFV charged lepton decays in the SM and SUSY GUTs and present the 
predictions for the five SUSY SO(IO) models considered. In Sec. IV their expectations for fi — e 
conversion are given. Sect. V concludes this paper. 

II. PREDICTIVE SUPERSYMMETRIC GRAND UNIFIED MODELS CONSIDERED 

We begin with a brief discussion of the general formalism on which the supersymmetric 50(10) 
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grand unified models are based. For all five models to be illustrated, the seesaw mechanism 
is of the conventional type I leading to normal hierarchies for the light Majorana neutrinos. The 
leptonic sector of the Yukawa superpotential at the GUT scale can then be written as 

Wy = N^Y^jLjHu + EfYf^LjHd + N^M^jN^, (4) 

where L represents the left-handed lepton doublet, A'^'^ the left-handed conjugate neutrino 
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singlet, and the left-handed conjugate charged lepton of one of the three 16 dimensional 
representations of 50(10). When the two Higgs doublets, and Hfi, acquire vacuum expectation 
values, the charged lepton and neutrino mass matrices are generated and can be written in 6 x 6 form 




(5) 



Mn Mr 

With the entries in Mr much larger than those in M^r, the light Majorana neutrino mass matrix 



is given by the well-known type I seesaw formula 



111 



= -MJjM^^Mn- (6) 

Note that the above Dirac 3x3 mass matrix entries. Me and M^v, are written in right - left 
order and appear in the 50(10) flavor basis. As such, the matrices can be diagonalized by the 
unitary transformations 

uI;^MeUel = diag(me, m^, m^), 
uIrMrUI,^ = diag(Mi, M2, M3), (7) 
Uj^M^UuL = diag(mi, m2, nis), 

by forming the Hermitian products, calculating the left and right unitary transformations of column 
eigenvectors, and phase rotating the complex symmetric Mr and matrices, so their mass 
eigenvalues are real. The PMNS neutrino mixing matrix [12,] is then given by Vrmns = U'^ejJJ^l- 
On the other hand for later use, it is convenient to transform to the bases where Me and Mr 
are diagonal and denoted by primed quantities. In this case 

Ml = -M'^M'-^M'^, 
KI^^Kl = diag(mi, ma, ms), 

where now Vpmns = U^^, since U'^^ is just the identity matrix. A comparison of the seesaw 
diagonalization matrices then reveals that the transformed Dirac neutrino matrix in the new 
basis is given by 

= uIj^MnUel, (9) 
in terms of the original matrix in the flavor bases. 
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With either basis, the left-handed neutrino PMNS mixing matrix can be written as Vpmns 
UpMNS^j where by convention [l^ 



UpMNS 



^ C12C13 S12C13 sise""^ ^ 

-S12C23 - Cl2S23Sl3e*'' C12C23 - Si2S23Sl3e*'^ S23C13 
\ S23S12 - Cl2C23Si3e*'^ -C12S23 - Si2C23Sl3e*'^ C23C13 / 



(10) 



in terms of the three mixing angles, 612, O23 and ^13; and the Dirac CP phase, d, in analogy with 
the quark mixing matrix. The Majorana phase matrix, 

$ = diag(e^^i, e'^\ 1), (11) 

written in terms of the two Majorana phases, xi X2j is required since an arbitrary phase 
transformation is not possible when one demands real diagonal neutrino mass entries in the trans- 
formation of M,y in Eq. dH) or in Eq. ([8]). 

With this background in mind, we now turn to a brief discussion of the five SO{10) models we 
have considered for our lepton flavor violation study. The 5*0(10) grand unification symmetry is 
an economical and attractive one [7], for all sixteen left-handed quark and lepton fields and their 
left-handed conjugates fit neatly into one 16 representation per family. Many models exist in the 
literature which differ from one another by their Higgs representation assignments and type of flavor 
symmetry imposed, if any. To appreciate this, it is of interest to note the following decompositions 
of the direct product of representations: 

16 16 = 10s e 120^0 126s, 

_ (12) 
16 16 = 10 45 210, 

where in the first product the 10 and 126 matrices are symmetric, while the 120 is antisymmetric 
under the interchange of family indices. 

Given this group structure for S'O(IO), there are two general classes of models which have been 



extensively studied. Those with Higgs in the 10, 126, 126, and possibly the 120 and/or 210, 
dimensional representations are often referred to as the minimal Higgs models [14] and lead to 
symmetric or antisymmetric matrix elements, or a superposition of the two. The advantage is that 
the couplings are renormalizable and preserve R-parity, but the latter representations are of rather 
high rank and disfavored in the string theory framework. The other class typically involves the 
lower rank Higgs representations such as 10, 16, 16,45, with some non-renormalizable effective 

n 

operators formed from them, [15] but R-parity is not conserved. They can lead to lopsided mass 
matrices for the charged lepton and down quark mass matrices, due to the SU{5) structure present 



in the 16 's. As such, one may anticipate that they will predict a higher level of lepton flavor 
violation than the first class. The SU{2)l triplet components in the Higgs representations in all 
five models considered have no electroweak VEV and hence lead to the conventional type I seesaw 
mechanism with the prediction of normal hierarchy for the light left-handed neutrino spectrum 
[l^ . We elaborate on each model in turn but give only the neutrino and charged lepton mass 
matrices. All successfully predict the observed quark structure and CKM mixings. 



A. Albright - Barr S'O(IO) Model 



This model, based on 5'O(10) with a U{1) x Z2 x Z2 flavor symmetry [171], is of the lopsided 
variety and has matrices of the following textures: 

^ / (5 5' ^ 

mu, Me = 



Mn 



( 1] 5n S'j^ 
6n -e 

where the parameters have the values rj = 1.1 x 10~^, 
0.00827e^ii9-4°, 6n - 



rriD, 



6 -e 
\6' a + e 1 J 
e = 0.147, a = 1.83, 6 = 0.00946, 6' 



(13) 



-1.0 X 10-^ 



-1.5 X 10" 



mu 



113 GeV, mo = 1 GeV. The 



right-handed Major ana mass matrix is given by 



/ c^r/^ —ber] arj 



Mr 



-ber] e — e 
\y arj —el 

where a = c = 0.5828i, b = 1.7670z, Ar = 2.35 x lO^'' GeV. One finds 

Ml ~ M2 ~ 4.45 X 10*^ GeV, M3 = 2.4 x 10^^ GeV, 
nil = 3.11 meV, 1712 = 9.48 meV, 772.3 = 49.13 meV, 

^21 



(14) 



Am^i = 8.0 X 10"^ eV^, Amj^ = 2.3 x lO^^ eV^, 



(15) 



sin2 26*23 = 0.99, sin^ 6*12 = 0.28, sin^ 6*13 = 0.0020. 



A value of tan/? = 5 assures that a ^ e and that the corresponding lopsided nature of the 23 
element of the down quark mass matrix gives a good fit to the results for the quark sector. The 
6^, elements of the Dirac neutrino matrix were added [18;] to the original model to give a 
better fit to baryogenesis arising from resonant leptogenesis involving the two lighter right-handed 
neutrinos [l^ . Evolution downward from the GUT scale to Mz has little effect on the mixing 
angles due to the small tan (3, the opposite CP-parity of the lighter two right-handed Majorana 
neutrinos, and their large hierarchy with the heaviest one; but serves to lower the values of the 
mj's and A'm'^2 ^^'^ ^^21 relative to their GUT scale values. 



B. Chen - Mahanthappa S'O(IO) Model 

This model, based on 50(10) with a SU{2) x Z2 x Z2 x Z2 flavor symmetry [2^, is of the 
minimal Higgs variety leading to symmetric entries in the mass matrices: 

/o a\ / ee-'^ 0^ 

vdsinP, Me = 



Mn 



be'" c 
a c 1 



vhcos /3, 



(16) 



ee 

ee*"^ -3/ 

ya c 1 J y I J 

where v = 174 GeV and the parameters have the values a = 0.00250, b = 0.00326, c = 0.0346, d = 
0.650, e = 0.004036, / = 0.0195, h = 0.06878, d = 0.74, cf) = -1.52. The solar and atmospheric 
neutrino mass squared differences, Am^^j = 8.14 x 10~^ eV^ and — 2.3 x 10~^ eV^ were 

used as input to determine the t = 0.344 and M3 = 6.97 x 10^^ parameters in the effective light 
left-handed Majorana neutrino mass matrix 

/ i ^ 



M, 



1 1 + 1" 
\t l+f 1 



(udsin I3y 
Mr 



(17) 



One then finds 



Ml = 1.09 X 10^ GeV, M2 = 4.53 x 10^ GeV, M3 = 6.97 x 10^^ Q^y, 
mi = 2.62 meV, m2 = 9.39 meV, = 49.2 meV. 
sin2 2l923 = 1-00, sin^ 6^12 = 0.27, sin^ 6^13 = 0.013. 



(18) 



For this model, a value of tan /? = 10 is used. The effect of evolution from the GUT scale to Mz 
is to raise sin^ ^12 and to lower sin^ 6*13 . 



C. Cai - Yu 50(10) Model 

I — I 

This model is based on 50(10) with an ^4 flavor symmetry [21]. The Higgs fields appear in 
six lO's, three 126's, three 126, and one 210 representations of 50(10), distinguished by their 54 
flavor assignments. Of the 14 pairs of Higgs doublets at the GUT scale, all but one pair are assumed 
to get superheavy. The charged lepton mass matrix is chosen to be diagonal, while the right-handed 
Majorana neutrino mass matrix is proportional to the identity matrix, so all three heavy neutrinos 
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are degenerate. The Dirac neutrino and charged lepton mass matrices are symmetric and given by 

^ ao — 2a2 — 3((io — 2^2) 05 04 ^ 

Mn = 05 oo + ai + 02 - 3(do + di + ^2) as 

Y a4 as ao - ai+ 02- 3{do - di + ^2) y 

^60-262-3(60-262) \ 

= 60 + ^1 + ^2 - 3(eo + ei + 62) 

^ 60 - 61 + 62 - 3(eo - ei + 62) y 



(19) 

where the parameters have the values ao = 18. 935+0. 000271681z, ai = -30.7989+0.0019887?, 02 = 
10.0361 + 0.00171701i, as = -2.99072 - 0.054757i, 04 = 0.554859 - 0.234705i, 05 = -0.066748 + 
0.008155i, 60 = 0.387756, h = -0.539649, 62 = 0.19327, do = 8.60218, di = -10.1912, ^2 = 
3.72519, eo = -0.0227734, ei = 0.0228717, 62 = -0.0115298, all in GeV. The common mass 
of the degenerate right-handed neutrinos is determined to be Mji = 2.4 x 10^^ GeV, so as to fit 
Arn^i = 2.6 x 10^'^ eV^ with the aid of the seesaw formula. The authors find 



(20) 



mi = 7.7 meV, m2 = 11.8 meV, ms = 50.7 meV, 
sin2 26'2s = 1.00, sin^ 6*12 = 0.29, sin^ 613 = 0.0029. 

A value of tan /? = 10 is used to evolve the quark and charged lepton masses, though the neutrino 
masses and mixings have not been evolved downward to the electroweak scale. 



D. Dermisek - Raby 50(10) Model 



22l |. The charged lepton and down 



This model is based on SO (10) with a family symmetry 
quark mass matrices have lopsided textures, but they are not so extreme as in the case of the 
Albright - Barr model: 



Mn 



e'uj 



( ° 

—e'uj 3euj 



-3e^cj 
-3e(T 
1 



\ 



vXsinP, Me 



^ e' -Se^cr^ 

—e' 3e —3ea 
^ 3eC 3e 1 j 



vXcosp, (21) 



where v = 174 GeV and the parameters have the values A = 0.64, e = 0.046, a = 0.83e' 



0.618j 



OMle^-^^^\ p = -0.053e"-^'*^*, e' = -0.0036, ^ = 0.12e^-^^^*, u = 2a/{2a - I). The right-handed 
Majorana mass matrix is diagonal with mass eigenvalues Mi = 1.1 x lO^'', M2 = —9.3 x 10^^, 
Ms = 5.8 X 10^^ GeV. These parameters were chosen by making use of the central experimental 



„0.767j 



values for Am^^ = 7.9 x 10"^ eV^ Amj^ = 2.3 x lO'^ eV^, sin^ (9i2 = 0.295 and sin^ 6*23 = 0.51 
at the time of writing. The authors find 

m^ = 3.7 meV, mo = 9.6 meV, = 49.2 meV, 

(22) 

sin^ (9i3 = 0.0024. 

The evolution has been carried out with tan P = 49.98 and ^0 = —6888.3 GeV down to Mz and 
then down to the 1 GeV scale for the light neutrino masses. 



E. Grimus - Kuhbock 50(10) Model 



This model is based on S'O(IO) with an Z2 flavor symmetry [23(]. The Higgs fields appear in one 



each of the 10, 120, and 126 dimensional representations of S'O(IO). The fermion mass matrices 
are generated by renormalizable Yukawa couplings of these Higgs fields to the three families placed 
in 16's. The Dirac neutrino and charged lepton mass matrices are written as linear combinations 



of the Yukawa couplings of the 10, 120, and 126 Higgs representations, respectively, while the 
right-handed Major ana matrix is proportional only to the latter: 

Me = H' + rLe''^i^G' -Se^'^^F' 

(23) 

Mr = rR^^F', 



G' 



\ 
/ 

V 
/ 

v 



10" 



-3 



(24) 



where the individual Higgs contributions to the mass matrices are given by 

'^ 0.716986 ^ 

H' = -40.6278 x 10~^ 

1114.41 J 

( -7.56737 \ 
7.56737 -36.8224 

36.8224 
-0.0966851 4.25282 ^ 
F' = 12.3136 X 10 

4.25282 -61.6491 

in the right-left order we have adopted. The coefficient parameters are taken to be rn = 
91.0759, rp = 297.758, r„ = 7.14572, = 1.33897, = 3008.88, = 2.90553 x 10^l^ Cd = 
19.66974°, Cn = -2.96594°, V'L = 6.24258°, = 179.85271°. A value of tan/3 = 10 is used to 



evolve the quark and charged lepton masses downward to the electroweak scale. One finds then 
that 

mi = 1.6 meV, mo = 9.2 meV, rn-'^ = 50.0 meV, 

(25) 

sm^ 26^23 = 1-00, sin^ 6*12 = 0.31, sin^ 613 = 0.00059. 

Although the results are impressive, the authors do note that 21 parameters have been introduced 
in order to obtain the results. Hence the model should be viewed as an existence proof that 
the lepton masses and mixings as well as the quark masses and mixings can be described in the 
framework of renormalizable couplings with only the three Higgs representations contributing to 
the Yukawa couplings. 

We close this section by summarizing the features of these models in Table I. Three of the models 



have similar predictions for sin^ 
Daya Bay reactor experiments 



0.0025, near the expected reach of the Double CHOOZ and 



24l ]. One of the models predicts a value near the present upper 



TABLE I: Higgs representations, flavor symmetries, and other noteworthy features of the five 50(10) SUSY 
GUT models considered in this work. 



Models 


Higgs Content 


Flavor Symmetry 


Mr (GeV) 


tan/3 


sin^ 6113 


Interesting Features 


AB 


10, 16, 16,45 


U{1) XZ2XZ2 


2.4 X lO^"' 

4.5 X 10* 
4.5 X 10* 


5 


0.0020 
(2.6°) 


Large Mr hierarchy with lightest 
two nearly degenerate leads to 
resonant leptogenesis. 


CM 


10, T26 


SU{2) X {Z2f 


7.0 X 10^2 
4.5 X 10^ 

1.1 X 10^ 


10 


0.013 

(6.5°) 


Large Mr hierarchy with heaviest 
more than 3 orders of magnitude 
below GUT scale; large sin^^ia. 


CY 


10, T26 


Si 


2.4 X 10^2 
2.4 X 10^2 
2.4 X 10^2 


10 


0.0029 

(3.1°) 


Degenerate Mr spectrum 4 orders 
of magnitude below GUT scale. 


DR 


10, 45 


Dz 


5.8 X 1013 
9.3 X 10" 
1.1 X 10^° 


50 


0.0024 
(2.8°) 


Mild Mr hierarchy almost 3 orders 
of magnitude below GUT scale. 


GK 


10, 120, 126 


Z2 


2.1 X 10^5 

4.2 X 10^'' 
6.7 X 10^2 


10 


0.00059 
(1.4°) 


Mild Mr hierarchy just 1 order of 
magnitude below GUT scale; 
rather small sin^ ^13 . 
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H- U^k Ul^ e- yj-. y^ Vk 

FIG. 1: Example of a Feynman diagram for fj, ^ ej with a neutrino mass insertion in the SM. 
limit placed by the CHOOZ experiment 



25|, while the fifth model predicts a value of 3 x 10~^ 
which is essentially beyond reach until a Neutrino Factory becomes a reality. We shall see what 
potential success charged lepton flavor violation experiments will have in further distinguishing the 
viable models. 

III. LEPTON FLAVOR VIOLATION IN RADIATIVE DECAYS 

We now turn to the subject of charged lepton flavor violation that can occur in the radiative 
decays, >e + 7, + and r — > e + 7. In the SM with the addition of three massive 

right-handed neutrinos, we observe that the individual lepton numbers, Le, and Lr, are not 
individually conserved. In radiative lepton decays, the flavor violation arises in one loop, where 
the neutrino insertion involves lepton flavor-changing Yukawa couplings of the left-handed and 
right-handed neutrinos; cf. Fig. 1. The branching ratio BR21 defined by the ratio of the rate for 
the /Lt — > 67 mode relative to the purely leptonic mode, /i v^j^eve-, is given by [26| 

2 

(26) 



BR21 = ^ 



~ 3a /" Am|^ '\ ^ ■ 2 ofl ^ in-54 
- T287 \-M^ ) 2(^12 ~ 10 , 

where the C/'s are elements of the PMNS mixing matrix. Hence in the SM, the expected branching 
ratio is immeasurably small, and the MEG experiment [27] looking for ^ — > 67 would be expected to 
give a null result. The present upper limit of 1.2 x 10~^^ was obtained by the MEGA collaboration 



28|. 



In SUSY GUT models on the other hand, the leading log approximations involve slepton- 

I i 

neutralino and sneutrino-chargino loops which contribute to the radiative lepton decays [6']; cf. 
Fig. 2. With more comparable heavy masses of the SUSY particles in the loops and lack of 
a Glashow-Iliopoulos-Maiani (GIM) mechanism, such a great suppression of the branching ratio 



does not occur. We shall work in the CMSSM 



29|, where the soft-breaking scalar and gaugino 



masses and trilinear scalar couplings are assumed to be universal at the GUT scale. The lepton 
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FIG. 2: Examples of Feynman diagrams for slepton - neutralino and sneutrino - chargino contributions to 
yu ^ 67 in SUSY models with slepton mass insertions. 

flavor violation then arises from evolution of the Yukawa couplings and soft-breaking parameters 

□ 

from the GUT scale down to the electroweak scale j3(|. 

Feynman diagrams for the LFV radiative decays in leading log approximation involve both 
neutralino - slepton {x^ — £) loops and chargino - sneutrino (x^ — 9) loops with the emitted 
photon attached to the internal charged slepton or chargino, respectively. Through evolution from 
the GUT scale, the LFV neutrino Yukawa couplings are induced primarily in the mass squared 
submatrix for the SU{2)l doublet sleptons, rn^iLL). We do not repeat the details here for this 
complicated calculation but rather refer the reader to the pioneering paper of Hisano, Moroi, Tobe, 
and Yamaguchi [s^ . We simply note that the radiative decay rate is given by 

r(^7 - C7) = ^rnl + P + \Af + A^^ \') , (27) 

where (n) and (c) refer to the neutralino and chargino loop contributions to the transition form 
factors Al and Aji connecting leptons of opposite chirality. The branching ratio for the flavor- 
violating decay mode relative to the flavor-conserving purely lepton mode is then 

BRii- - 17 = ^(I^Lp + \Ar\'). (28) 

In the leading log approximation with the largest contribution coming from the left-handed 
slepton mass matrix, the branching ratio is given by 

^Rji = T^li^L^tan^P, (29) 

where 

«Lh = -^rnl{3 + Al/ml)Yl,log n,. (30) 

12 



with the Yukawa couplings specified in the lepton flavor basis and the right-handed Majorana 
matrix diagonal, so is just the kth heavy right-handed neutrino mass, while Mq is the GUT scale 
typically equal to 2 x 10^^ GeV and rris is some typical SUSY scalar mass. Petcov and collaborators 



3l[ have shown that the full evolution effects as first calculated in [30] can be extremely well 



approximated by Eq. ([22D , if one sets 

ml ~ Q.f>mlMli^{ml + Q.^Ml/^f. (31) 

While the branching ratio for ^ e+7 is well approximated by the above equations, the branching 
ratios for r ^ /i + 7 and r — > e + 7 must be scaled by the branching ratios for the flavor-conserving 
leptonic modes relative to their total decay rates. 



We see that the MEG experiment [27[ has no chance of observing a positive signal in the 
SM for the ;U — > 67 decay channel; however, the situation will be totally different for the SUSY 
GUT models we have chosen to consider. In the CMSSM with universal soft-breaking parameters 
mo, and for a given tan/3 and sgn{fi), we consider the following correlations between 

the branching ratio and the soft-SUSY breaking parameters: 

(a) the branching ratio vs. Mi ^2 fo^' fixed Aq = 0, for example, with different choices of mo; 

(b) the allowed parameter space for Ao/mo vs. M1/2, for specific branching ratio ranges; 

(c) the branching ratio for r — > /i7 vs. that for ^ — >■ ej, on a log-log plot, which are related by 

logBR(r ^ n-f) = logBR(/i ^ ej) + log(BR32/BR21) + logBR(r ^ i^rfJ-j^f,) 



log BR(/i 67) + log 



{YjLY^)2i 



0.757, 



where L = log(M||/M^), and again the Yukawa matrices are converted to the lepton flavor 
basis with Mr diagonal. Due to the factorization of the soft-breaking parameters and the 
GUT model parameters in the approximate Eq. (j29p . the slope is unity, and the intercept is 
just the sum of the last two terms after correcting for the r z^r/UP^ branching ratio. The 
length of the straight line segment depends on the range of the soft-breaking parameters 
chosen. A similar plot can be made for the branching ratio for r — > 67 vs. that for fi — > ej. 



13 



For all three types of plots we have imposed the following soft parameter constraints [13l |: 



For tan/? = 5, 10 : 


niQ : 


50 - 


400 GeV 




Ml/2 ■■ 


200 - 


1000 GeV 




Ao: 


-4000 - 


4000 GeV 


For tan /? = 50 : 


mo : 


500 - 


4000 GeV 




Ml/2 ■■ 


200 - 


1500 GeV 




Ao : 


-50 - 


50 TeV 



(33) 



In addition it is desirable to impose WMAP dark matter constraints in the neutralino, stau or 
stop coannihilation regions 



33l |. where the lightest neutralino is the LSP. These more restrictive 



constraints are wel 
soft gaugino mass 



described by the quadratic polynomial for the soft scalar mass in terms of the 



mo = co + ciMi/2 + C2Mf^^, 



(34) 



Ci 



Ci{Ao, tan/3, sgn{fx)). 

where . bounded since M,, is bounded. If M,,, is too s.all, the ptesen. experimental bound 

on the Higgs mass M of mfi > 114.4 GeV may be violated or the neutralino rehc density in the 
early universe will be too small, while if M1/2 is too large the neutralino relic density will be too 
large. We shall impose both limits on the DM constraint conditions for various values of Ao and 
find that both lower and upper limits are placed on the branching ratios for each model. 
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FIG. 3: Branching ratio predictions for /i ^ e + 7 in the Albright - Barr model -with tan/3 = 5. In (a) Aq 
is set equal to zero, "while in (b) all three parameters, mg, M1/2, Aq, are allo"wed to vary. 
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FIG. 4: Branching ratio predictions for /i — > e + 7 in the Chen - Mahanthappa model with tan [3 = 10. 



For each model in Figs. 3(a), 4(a), 5(a), 6(a), and 7(a), we plot curves with constant values of 
mo indicated for the BR21 branching ratios as functions of Mi 12 with the trilinear coupling = 
at the GUT scale. The curve which cuts across the rriQ curves in each left-hand plot represents 
the WMAP dark matter constraint for this value of Aq. The horizontal broken line indicates the 
present experimental upper limit for this branching ratio 

In Figs. 3(b), 4(b), 5(b), 6(b), and 7(b), we allow Aq to depart from zero and show scatterplots 
of Aq/itio vs. Ml/2- The points are color-coded as indicated according to the branching ratio 
intervals in which they fall, with all points below the present upper limit on the ^ 67 branching 
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FIG. 5: Branching ratio predictions for fi 



7 in the Cai - Yu model with tan /3 = 10. 
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FIG. 6: Branching ratio predictions for — > e + 7 in the Dermisek - Raby model with tan /3 = 50. 

ratio. The soft parameter constraints in Eq. (j33p have been imposed for the Monte Carlo selection 
of points. The continuous curves represent the WMAP dark matter constraints and are drawn in 
steps of 500 GeV from ^0 = -2.5 TeV to Aq = 2.5 TeV. 

It is clear that the predicted branching ratio decreases as the universal soft gaugino mass at 
the GUT scale increases along the constant scalar mass curve. The Grimus-Kuhbock (GK) and 
Albright-Barr (AB) models will be probed first and then the Dermisek-Raby (DR) model by the 
MEG experiment, while the other two models are essentially beyond reach, if = as depicted. 



(a) (b) 




FIG. 7: Branching ratio predictions for /i ^ e + 7 in the Grimus - Kiibock model with tan/3 = 10. 
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One sees that higher values of BR21 are predicted for a given value of as |74o/mo| increases. 
For the AB and GK models the experimental branching ratio greatly limits the allowed ranges of 
Ao/mo, while for the Chen-Mahanthappa (CM), Cai-Yu (CY), and DR models the dark matter 
constraints limit the allowed ranges of Ao/rriQ. In any case, the minimum predicted BR21 branching 
ratio occurs for Aq = 0. 

One can also present similar scatter plots for the r — > ^7 and r — > 67 decay modes. Since 
all three decay modes are intimately related in each model through the corresponding logarithmic 
terms such as that in Eq. (j32p . the same scatter points will appear with only the color-coding 
changed (assuming one imposes the BR21 experimental limit for each plot). 

Instead, we present two log-log plots for the BR32 and BR31 branching ratios against that for 
BR21 in Figs. 8 and 9, where ^0 = has again been imposed. The thin line segments for each 
model observe the soft parameters constraints imposed, while the heavier line segments observe 
the more restrictive WMAP dark matter constraints. The vertical dashed line reflects the present 
BR21 bound, while the horizontal dashed line refers to the present BR32 or BR31 experimental 
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FIG. 8: Branching ratio predictions for r —>/_* + 7 vs. branching ratio predictions for /i ^ e + 7 in the five 
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more restrictive WMAP dark matter constraints apply for the thick line segments. The present experimental 
constraints are indicated by the dashed lines. 



17 



10 
10 
10 
10 

+ 10" 
^ 10 
^ 10 



-06 



-07 



-08 



-09 



10 



12 



■13 



10 



10 



10 



14 



■15 



16 



10' 



10 



I I I mill — I I 1 1 Mill — I I I mill — I I 1 1 iiiij — I I 1 1 iiiij — I I I iiiii| — I I 1 1 iiiij — I I 1 1 nil 



— AB model 
. CM model 

— CY model 
DR model 

— GK model 




I I I I I I I I 



■18 

10-^' 10-'^ 10-'^ 10-'" 10-'^ 10-'^ 10-'^^ 10-'^ 10-" 10-'° 10-°' 10-' 

BR( \i-> e + Y ) 



FIG. 9: Branching ratio predictions for t — )■ e + 7 vs. branching ratio predictions for /i ^ e + 7 in the five 
models considered. The soft SUSY breaking constraints imposed apply for the thin line segments, while the 
more restrictive WMAP dark matter constraints apply for the thick line segments. 



limit, respectively [35|. It is clear from these two plots that the ongoing MEG experiment stands 
the best chance of confirming the predictions for or eliminating the GK and AB models. Even with 
a super-B factory 



36( 1 ■ the present experimental bounds on the BR32 and BR31 branching ratios 
can only be lowered by one or two orders of magnitude at most. 

But recall that the line segments apply for the special case of Aq = 0. If one allows Aq to depart 



TABLE 11: Ratios of the branching ratios for the lepton flavor violating r decays and fi — e conversion on 
Ti relative to that for fi decay. 



Models 


BR(t - 


M7)/BR(/i ej) 


BR(t - 


e7)/BR(/x 67) 


BR(^ + Ti- 


e + Ti)/BR(/i ^ 67) 


AB 




16.7 




0.09 




0.33 


CM 




1.3 X 10^ 




171 




0.11 


CY 




400 




6.5 




0.11 


DR 




3.3 X 103 




61.0 




0.026 


GK 




10.0 




1.0 




0.12 
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from zero, the line segments will slide diagonally upward and toward the right along their presently 
depicted positions by amounts that can be estimated from Figs. 3-7. Hence only the lower limits 
on the branching ratios are robust in Figs. 8, 9 and 11. However, it is clear that ratios of the 
branching ratios remain fixed for each model for any allowed Ao/mQ. We present these ratios for 
the T and /x — e conversion branching ratios relative to the fi ^ branching ratio in Table II. The 
spread in numbers for the ratios in different models appears to be greater than that anticipated 
by the authors of Ref. [s^] for the class of models considered here. 

In Table HI we summarize the relevant findings from our study of the five models. The branching 
ratio ranges apply for the Aq = case and with the stricter WMAP dark matter constraints 



TABLE III: Summary of the relevant results for the five SO{10) SUSY GUT models considered in this work. 
The present experimental upper limits for the branching ratios are indicated in the second line, while the 
third line of the table gives the projected upper limit reaches for the Meg experiment, Super-B factory, and 
next generation /i — e conversion experiment. 



Models 
Expt. 
Limits 


sm &i3 


Mr S 

(GeV) 


tan p 


\ A / 1 
|^o/»7to|max 


"r>T> o 1 / \ 

BK21(/i — ^ 67) 
< 1.2 X 10-" 
^ < 10-" 


BK32(t — > ii-f) 
< 4.5 X 10-8 
^ < 10-9 


BH(/i + ii-—i-e + ii) 
< 4 X 10-12 
^ < 10-" 


AB 


0.0020 
(2.6°) 


2.4 X lO^"^ 

4.5 X 10^ 
4.5 X 10^ 


5 


5 


(0.2-9) X 10-12 


(0.03-1) X 10-" 


(0.03-2) X 10-12 


CM 


0.013 
(6.5°) 


7.0 X 10^2 
4.5 X 10^ 

1.1 X lO'' 


10 


12 


(0.02 - 4) X 10-" 


(0.02-5) X 10-" 


(0.01 - 3) X 10-" 


CY 


0.0029 
(3.1°) 


2.4 X 10^2 
2.4 X 10^2 
2.4 X 10^2 


10 


12 


(0.02-5) X 10-" 


(0.04-9) X 10-" 


(0.03-6) X 10-" 


DR 


0.0024 
(2.8°) 


5.8 X 10^3 
9.3 X 10" 
1.1 X 10^° 


50 


2.5 


(0.05-8) X 10-" 


(0.02 - 3) X 10-9 


(0.01 - 2) X lO-i-i 


GK 


0.00059 
(1.4°) 


2.1 X 10^5 

4.2 X 10^^ 
6.7 X 10^2 


10 


2 


(0.4-80) X 10-" 


(0.004- 1) X 10-8 


(0.02-5) X 10-11 
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imposed. It is clear that the five predictive 50(10) SUSY GUT models considered have very 
representative right-handed neutrino mass spectra and predictions for sin^^ia. The CM, DR, and 
GK models have massive hierarchical spectra with M3 ranging from 10^'^ to 10^^ GeV. The CY 
model, on the other hand, has a degenerate spectrum with Mr ~ 3 x 10^^ GeV, while the AB 
model has degenerate Mi and M2 which can lead to resonant leptogenesis. The CM model has 
a relatively large sin^ ^13 prediction which should be observable at the upcoming reactor neutrino 
experiments. Double CHOOZ and Daya Bay. The AB, CY and DR models have similar predictions 
for sin^ ^13 which will make observation of Ve oscillation somewhat marginal at those reactors 



3a] without a SuperBeam source. 



and in the proposed NOz/A and T2K long-baseline experiments 
For the GK model, the observation of such a low sin^ ^13 prediction would only take place with 
a Neutrino Factory. But it is clear from Table III and the previous figures that the GK and AB 
models will be tested first with the MEG experiment. From our discussion it is then clear that the 
LFV branching ratios are more sensitive to large M3 than to large ^13 in the models considered. 
In previous generic studies of SUSY GUT models, the rare branching ratios were nearly equally 
sensitive to each of the two parameters 39l |. 



IV. LEPTON FLAVOR VIOLATION IN ^ e CONVERSION 

It is also of interest to consider lepton flavor violation in the /i to e conversion process in 
Titanium, /x + Ti ^ e + Ti. While there is no such ongoing experiment, preliminary discussion is 
underway to propose one which will lower the present limit. We shall see that the conversion rate 
predictions relative to the muon capture process, ^ + Ti ^ + Sc, are such that all five models 
considered in this paper can potentially be eliminated, if no signal is observed. 

While predictions for the conversion process in the SM are infinitesimally small, large enhance- 
ments in the CMSSM framework again occur by virtue of massive super partners appearing in loop 
diagrams involving 7, Z and Higgs penguins and boxes. Of these, Arganda, Herrero and Teixeira 
have shown that the 7 penguin contributions dominate the others by at least two orders of magni- 
tude 3] • The slepton-neutralino and sneutrino-chargino loop contributions to the 7 penguins are 
shown in Fig. 10. The virtual massive A^'^ and N'^ with their Yukawa couplings again appear in 
the slepton loops along with a higgsino or Higgs particle, respectively. 

The complex formulas for legton flavor violation in this process were first derived by Hisano, 



Moroi, Tobe, and Yamaguchi 



30|. If we restrict our attention to the 7 penguin contribution as 
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FIG. 10: Examples of Feynman diagrams for slepton - neutralino and sneutrino - chargino contributions to 
— e conversion in SUSY models with slepton mass insertions. 



suggested by Arganda et al. [4d[, one finds the /i — e conversion rate is given by 

r(^ - e) = ^oC^Zt^Zm^Fm^ {\A\ - A^f + |Af - A\f\ , (35) 

where here A^'^ are form factors for the vertices connecting leptons of equal chirality, while A2'^ are 
combinations of the electric dipole and magnetic dipole transition form factors previously denoted 



by ^L,_R- For a 22^^ target, Z^ff = 17.6 and the nuclear form factor is F{q'^ ~ — 0-54 30(]. 



In the case of the conversion process, we have explicitly carried out the full evolution running 
from the GUT scale to the Z scale. The — e conversion branching ratio is then obtained from 
the conversion rate above by scaling it with the /i capture rate on Ti, which is quoted in 4l| as 
(2.590 lb 0.012) X 10^ sec~^ with the present experimental limit on the conversion branching ratio 
found to be < 4 X lO"^^ 

In Fig. 11. we show a plot of the fi — e conversion branching ratio vs. the /i — > 67 branching 
ratio for each of the five models considered. We have limited the line segments by applying the 
WMAP dark matter constraints of Sec. III. It is clear that the GK and AB models would be 
tested first, followed by the DR, CY and CM models. In fact, a first generation fj, — e conversion 
experiment may be able to reach a branching ratio of 10~^^, while a second generation experiment 
may lower the limit from the present value down to 10^^^ {42^. If such proves to be the case and 
no signal is seen, all five models will be eliminated. Hence the conversion experiment is inherently 
more powerful than the MEG experiment looking for /u — > 67 which is designed to reach a level of 
10^^'^ — 10^^^, sufficient only to eliminate the GK and AB models. The caveat, of course, is that 
MEG is now starting to take data, while no new conversion experiment has been approved to date. 
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FIG. 11: Branching ratio predictions for fi — e conversion vs. branching ratio predictions for /i ^ e + 7 
in the five models considered. Only thick line segments are shown reflecting the apphcation of the more 
restrictive WMAP dark matter constraints. Note that the predictions for the CM and CY models overlap. 



V. CONCLUSIONS 



There have been many theoretical models constructed which aim to explain the neutrino masses 
and mixing patterns. While the predictions for the value of ^13 in these models may provide a 
way to distinguish some of these models, as shown in the survey over 63 models we have carried 
out in Ref. [5|], the rare LFV processes, such as /i — > e + 7 and fi — e conversion can provide an 
even more sensitive way to disentangle these models. We have investigated these rare processes 
in five SUSY SO(IO) models that are currently still viable and highly predictive, making use of 
the allowed parameter space for the soft SUSY breaking parameters in the CMSSM framework. 
Utilizing the WMAP dark matter constraints, lower bounds on the branching ratios of these rare 
processes can be placed, and we find that at least three of the five models considered give rise 
to prediction for ^ e + 7 that will be tested at MEG. More interestingly, the next generation 
/i — e conversion experiment should be sensitive to the predictions of all five models, making it an 
even more robust way to test these models. While generic studies have emphasized the important 
dependence of the branching ratios on the reactor neutrino angle, ^13, and the mass of the heaviest 
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right-handed neutrino, we find the latter to be by far the more significant in the models tested. 
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43] informed us that a suppression 



Note added.-Aftei completion of this work, the authors of Ref. 
factor should be applied to the lepton flavor violating branching ratios calculated in our paper. This 
suppression arises from a QED radiative correction to the effective dipole operators responsible for 
the rare decays. The multiplicative factor in question is given by 1—S, where 6 = (So/tt) log(A/mj), 
A ~ 250 — 1000 GeV is the sparticle mass scale responsible for the LFV, and mj = mfj_ or ttIt- for 
the rare radiative muon or tau decays. For ^ — > 67, the 5 correction amouts to roughly 0.15, so 
the branching ratio is about 0.85 times the rate given in the paper. The authors thank Andrzej 
Czarnecki and Ernest Jankowski for pointing out their earlier work to them. 
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